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We investigate spin transport in InSb/InAlSb heterostructure using the Monte Carlo 
approach, generalized by including density matrix description of spin for taking spin 
dynamics into account. In addition to the dominant Dyakonov-Perel (DP) mechanism 
for spin control and relaxation, we consider magnetic interaction between electrons 
which assumes importance due to high electronic Lande g-factor in the material. It is 
found that while the effect of magnetic interaction is not important at low densities, 
it reduces the spin relaxation length by as much as 50% at higher densities. We also 
present a calculation which elucidates the suppression of decoherence attributed to 
wave vector dependent magnetic field in the DP relaxation mechanism. We note that 
magnetic interaction is a general relaxation mechanism which may assume importance 
in materials with high electronic Lande g-factor. 
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I. INTRODUCTION 



The electron spin dynamics has been studied extensively due to the prospects spin based 
devices hold for future. In addition to offering improvements over the contemporary de- 
vices, Spintronics promises devices with entirely new functionalities^ - -. While metal based 
Spintronics has grown into a somewhat mature field over the past few years owing to several 
succesful applications of magnetic tunnel junctions, semiconductors still remain to make an 
impression. The importance of integrating spintronics with contemporary semiconductor 
technology cannot be overstated. 

Various III-V compounds are being investigated for the material to be used in semicon- 
ductor spintronics. The narrow gap semiconductors InSb and InAs offer distinct advantages 
over their wide gap counterparts 4 . Spin transport in both these semiconductors has been 
studied to find InSb better of the twc£. InSb has the lowest effective mass amongst III-V 
semiconductors, high mobility, and a stronger spin-orbit coupling. Further, the Rashba 6 
term has been found to dominate the Dresselhaus-^ term in the spin-orbit interaction at 
high densities 4 § . This, in principle, must offer a better control over the spin dynamics in 
InSb as Rashba coupling can be tuned by electric field. These factors make InSb a very 
promising candidate as semiconductor material for Spintronic applications in general, and 
for Datta-Das transistor- in particular. 

Of the three spin relaxation mechanisms in semiconductor systems, Bir-Aronov-Pikus 
mechanism^ is relevant for p type semiconductors only and Elliot- Yafet mechanism^ 1 - is 
dominated by Dyakonov-Perel (DP) mechanism^ in high mobility samples^ 1 ^. We propose 
yet another pathway to spin relaxation that gains importance in semiconductors with high 
g-factor. The electronic magnetic moment gives rise to an actual magnetic field which plays 
an important role in addition to the effective magnetic field attributed to spin-orbit coupling. 
This effect is proportional to g 2 and is three orders of magnitude higher in InSb (g ~ —50) 
as compared to other materials of interest such as GaAs (g ~ — 0.44)^. The motivation 
for considering magnetic interaction is derived from the observation that in ferromagnets, 
magnetic interaction leads to a phenomenon as important as the formation of domain walls. 
Hence, one is inclined to believe that magnetic interaction ought to play some role in restoring 
the spin space equilibrium in non magnetic systems for which g 2 becomes comparable to 
relative permeability in typical ferromagnets. 
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In this paper, we study the effect of electron-electron magnetic interaction on spin trans- 
port in high g-factor semiconductors using Monte Carlo simulations^^'i 6 -. Owing to the 
reasons discussed above, we study InSb/InAlSb heterostructure as a model system. How- 
ever, any given system can be simulated using our routine. The spin dynamics in the system, 
owing to its high mobility, is dominated by DP mechanism. We calculate the magnetic field 
attributed to non equilibrium spin concentration in a mean field approximation by assuming 
uniform distribution of charge and spin along the width of the channel. This magnetic field 
adds up to the pseudo field attributed to spin-orbit coupling to produce the net field ex- 
perienced by electrons. We have obtained dependence of spin relaxation length on electron 
density. A similar behavior has been observed in experimental studies on InSb 17 . In the 
section on results, we present a calculation which throws some light on the need to include 
magnetic interaction effect in addition to DP mechanism. Finally, we conclude with a dis- 
cussion on the requirement of a more accurate calculation than the mean field approach 
considered in this work. 



II. MODEL 

A detailed account of the Monte Carlo method used for studying spin dynamics in III-V 
semiconductors has been given elsewhere^^J^. We shall only enlist the essential features of 
the simulation. The dynamics has been assumed to be confined to the lowest subband. We 
treat the band structure via the nonparabolic approximation due to Conwell and Vassel^ 1 ^. 

E(l + aE) 

Here, a is given by the following expression^. 

Eg is the bandgap, m the effective mass of electron and is the mass of free electron. 
The scattering mechanisms incorporated in the simulation are acoustic phonon scattering, 
optical phonon scattering and ionized impurity scattering. The ionized impurity scatteing 
calculation considers the Conwell and Weisskopf approach 16,20 . We now discuss the boundary 
condition used in the system. In a realistic situation, there is a constant flux of electrons, 
independent of the applied voltage, from both electrodes. At zero voltage and steady state, 
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the inward flux from electrodes is exactly balanced by the outward flux from the channel. We 
exploit the fact that the outward flux from the channel changes marginally on application 
of a reasonable voltage. Hence, we choose to reintroduce the electrons, that fly into source, 
at the drain end and vice versa. This boundary condition distorts the spin polarization 
profile from exponential decay for a small channel device and presents a realistic picture, 
which allows us to take into account effects like spin selective transmission of electrons across 
drain boundary. 

Source Channel Drain Source Channel Drain 
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FIG. 1: Depiction of electronic flux in steady state at the two electrodes, Source and 
Drain, at (a) zero voltage and (b) small positive voltage Vds- Case (c) depicts the 
approximation that we consider in our calculations. Since 5<j) is very small, cases (b) and 

(c) are nearly the same. 

It is not possible to compute the complete 2D magnetic field intensity variation across 
the device for a given distribution of electrons at each simulation step. Hence, we adopt a 
mean field kind of approach towards its calculation and ignore any variation in magnetic field 
along the width of the device. For the purpose of simulation, the device is divided into a grid 
along the length of the channel. The magnetic field is calculated in two parts - interaction 
within the grid element and interaction between different grid elements. Interaction within 
the grid element is calculated assuming all the electrons in it to be spaced uniformly along 
width of the device and the net magnetic moment to be distributed equally amongst them. 

rtintragrid fin 1 N ^° ^x,y,zN 1 ( „s 

B *,y,z = (-lx,2„,-l z ) — m 2 2_,rf ^ 
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Coordinate system has been chosen such that x and y are along length and width of the 
device respectively. Here (—1,2,-1) corresponds to multiplicative factors for x,y and z 
components respectively N denotes the total number of electrons in the grid element, 
M x<yjZ denotes respective components of total magnetic moment due to all the electrons in 
the grid element and W is the width of the device. 

For calculation of inter grid element interaction, magnetic field due to each grid element, 
at other grid elements, is calculated assuming the net magnetic moment to be distributed 
uniformly along width of the device. It is a simple task of evaluating some integrals to find 
the required magnetic field intensity. Ignoring edge effects, magnetic field intensity due to 
a grid element 1 at a grid element 2 is given by the following expression. 

Tjinterqrid /i i \ A*0 M xz 1 

*•* = (1 " ~ lz) 27 ~W ^ (4) 
Symbols have their usual meanings (See Eq. (j3J)) and 7*12 is the distance between the two grid 
elements. It is observed that y component of intergrid interaction vanishes. The expressions 
above enable us to calculate magnetic field intensity along the channel length for a given 
spin polarization profile, which combined with spin-orbit coupling dictates spin dynamics of 
the electrons. 

In the range of densities considered in this work, Rashba term dominates spin-orbit 
coupling 4 . Hence, a major contribution to spin dynamics comes from the following Rashba 
Hamiltonian. 

H R = a (k y a x - k x a y ) (5) 

We consider Dresselhaus term in its linear approximation, the Hamiltonian for which is as 
given below. 

H D = (3 < k 2 z > {kyOy - k x a x ) (6) 

Addition of the Zeemann term takes care of the magnetic field intensity in a semi classical 
way. This method is approximate as a complete solution requires consideration of vector 
potential in the coordinate space Hamiltonian, which leads to further complications as free 
particle wave is no longer the solution. However the condition that magnetic field intensity 
(B) is small enough so that the momentum relaxation time is much smaller than m/eB en- 
sures that the effect of magnetic field on the coordinate space wavefunction remains minute. 
We deal with the effect of magnetic field in coordinate space by adding a Lorentz force term 
to coordinate space dynamics. 



III. RESULTS AND DISCUSSION 



The structure probed in simulations is similar to the one experimentally studied by 
Gilbertson and co workers 4 . The thickness of InSb layer in the considered device struc- 
ture is taken to be 20nm while both length and width are taken as 300nm. The choice of 
longitudinal dimensions is influenced by the inability to simulate more than 2 x 10 5 elec- 
trons. The values of a and < k 2 z > were taken as 1.5 x lO~ n eVm and 3 x 10~ 12 eVm 
respectively 4 -. The dynamics is assumed as being confined to the lowest subband, the energy 
for which was calculated as 0.1 eV. The various material and scattering parameters were 
adopted from a standard manual on Monte Carlo simulations- 1 . A remote doping scheme 
such as an adjacent doped layer allows for higher mobility and hence an effective impurity 
density of 2 x 10 n cm~ 2 is assumed to obtain a reasonable value of mobility. A basic time 
step of 0.1 fs was decided and electrons were evolved for 2 x 10 5 steps to achieve a steady 
state. Data for time averaging was recorded for the last 10000 steps only. The validity of 
simulation was thoroughly tested by reproducing already reported results^. 
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FIG. 2: \S\ along channel length. Thick lines denote the case when magnetic interaction is 

not accounted for in calculations, (a) Legend indicates electron density in cm -2 and 
electrons are polarized along z direction at the time of injection, (b) Legend indicates the 

polarization direction of electrons at the time of injection and n = 16 x 10 13 cm -2 . 
Maximum reduction in spin relaxation length of more than 50% is observed for the case of 

initial polarization along x. 
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Spin transport was simulated assuming 100% injection efficiency with Vds = 0.05V for 
different electron densities. Room temperature was assumed in all simulations. The results 
are reported in Fig. (j2aj) . We note that spin relaxation profile is independent of electron 
density when DP mechanism alone is considered. However, magnetic interaction is strongly 
dependent on density with its effect becoming reasonable at a density of 4 x 10 13 cm~ 2 . The 
decrease in spin relaxation length due to magnetic interaction is almost 30% for a density of 
16 x 10 13 cm -2 . It is to be noted that data close to drain was ignored while fitting in order 
to nullify the boundary effect. 

Anisotropic nature of the spin-orbit Hamiltonian and magnetic interaction leads to dif- 
ferent spin dynamics for different initial polarization of injected electrons. We compare the 



spin dynamics, with and without taking magnetic interaction into account, in Fig. f!2bj) . 
Our simulations show that the effect of magnetic interaction is most pronounced for the 
case when electrons are initially polarized in x direction. We compare the complete spin 
dynamics for this case in Fig. ( l3al) . We note that spin profile along the channel is modified 
substantially by the inclusion of magnetic interaction in our calculations. 
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FIG. 3: Initial polarization of electrons is along x and n = 16 x 10 13 cm~ 2 . (a) S x , y ,z along 
the channel length. Thick lines denote the case when magnetic interaction is not accounted 
for in calculations, (b) Magnetic field intensity variation along the channel length. 



We now discuss a calculation which brings out a somewhat diminished role of k depen- 
dent effective magnetic field in causing spin decoherence. The calculation also points to a 
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reduced decoherence attributed to variability of electron trajectories. We show that spin- 
orbit coupling is more important as a control mechanism than as a decoherence mechanism, 
and hence the need to consider other relaxation mechanisms for accurate spin dynamics. 
This result is in our favor as control via spin-orbit coupling is at the heart of many spin- 
tronic devices^. In what follows, we ignore the weak Dresselhaus term (Eq. (151)) at first 
and comment on its significance later. Considering the spin-orbit coupling alone, the spin 
dynamics is given by the following equations. 

da x 2a 



-r-k x cr z (7) 

2a , /nN 
- — k y a z (8) 



dt h 
dcjy 2a 

dt h 
da z 2a 

= -j-(k y a y + k x <j x ) (9) 

In the equations above, it is deemed as understood that u x . ytZ (which stand for Pauli matrices) 
and k XtVtZ represent their respective expectation values. Let us focus our attention on Eq. 
(I7j). Writing k x in terms of v x , integrating by parts and substituting da z /dt from Eq. ([9]), 



2am 
o x = o x - — 



2am , 

a z x - I (VyCTy + v x a x ) xdt 



(10) 



h 2 Jo 

In writing the above expression, it was assumed that at t — electron was introduced at 
x = with its initial polarization in the x-z plane. cr° and a z represent respective initial 
values. Now, integrating Eq. and substituting for <j z in the above equation, 



q 2amxa z Aa 2 m 2 
o x = a x - 2 



x 



t ft 



'vyOy + v x a x ) dt- x {vyOy + v x a x ) dt 
o Jo 



We pause here to build the background for discussion on the significance of expression derived 
above. Firstly, we note from simulation results (Fig. ( |3aj) ) that there is little change in o y 
from its initial value. This is a consequence of a strong electric field along the x direction 
which results in > \k y \ on an average, which causes the effective magnetic field due 
to spin-orbit interaction to be almost aligned along y direction. Hence, we are justified in 
treating a y as almost constant. Since a 2 . + a 2 + a 2 z = 1 and with a y fixed, only one of a x 
and a z is independent. Now we are ready to discuss Eq. (11 II) . We note that the first two 
terms on the right hand side are same for all electrons irrespective of their velocities or the 
paths they take to reach a particular x. Hence, the k dependent Hamiltonian does not lead 
to decoherence to the first order in spin-orbit coupling parameter a. The first order term 
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provides us with the much needed control of spin dynamics through spin-orbit interaction. 
The second order term, which is relatively small, is the one responsible for spin decoherence. 

Physically, this can be understood with the help of a simplified example. Consider an 
electron that moves from a point A to another point B and back to point A (Fig. flU)). As 
the electron moves from A to B, it undergoes precession about an effective magnetic field 
and thus suffers a change in its polarization. Now as it reverses the direction of its motion 
to move back to point A, the direction of effective magnetic field is reversed as well and 
the electron begins to precess in the opposite direction. The faster it moves, the faster it 
precesses so that when it comes back to point A, it has undergone almost the same change 
in polarization in the reverse direction. On the other hand, magnetic interaction does not 
have any coherent effect as the precession angle depends upon the path traversed by the 
electron, which is different for all of them. Further, the magnetic field is stronger in the 
region with higher spin polarization. 




A 



FIG. 4: Physical demonstration of coherent first order term. 

It should be noted that inclusion of Dresselhaus term in this calculation yields an addi- 
tional term (3 < k 2 z > y added to ax in the first order coherent term of Eq. ( TTT1) . The depen- 
dence of polarization on y coordinate can lead to substantial decoherence after summing up 
over all the electrons at a given x. This result elucidates the advantage in suppressing the 
importance of Dresselhaus term for better control of spin dynamics in the device. Hence, 
using a high electron density in the device might be useful for better spin control 4 in addition 
to increasing the speed. 

We conclude this section by emphasizing the need for a better 2D model than the mean 
field model considered for magnetic interaction in this work. In our calculations, all the elec- 
trons in a grid element experience the same magnetic field irrespective of their y coordinate 
and their distribution within the grid element. This is certainly not true and leads to an 
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underestimation of the actual decoherence rate. As a result, in our simulations, the density 
at which magnetic interaction effects gain importance comes out to be somewhat high. In 
actual devices, magnetic effects in InSb may set in at much lower electron densities. The re- 
ported magnetic effects can be responsible for the observed dependence of spin decoherence 
rate on electron density in InSb^. 

IV. CONCLUSION 

We have discussed a new pathway to spin decoherence in high ^-factor semiconductors. 
Our simulations show an increased spin decoherence attributed to electron-electron magnetic 
interaction. We have carried out an approximate calculation which elucidates the weakness 
of DP spin relaxation mechanism in the context of Datta-Das transistor. It also manifests 
the advantage of a diminished importance of Dresselhaus term in the overall spin-orbit 
interaction. InSb appears to be an appealing candidate for semiconductor based spintronics 
due to low electron effective mass and relatively strong spin-orbit coupling. At the same 
time, a high Lande ^-factor brings about additional effects. The electron-electron magnetic 
interaction leads to a shorter spin relaxation length. Hence, InSb might not be a suitable 
material for spin transport devices. 
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